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, ' In this paper, the block processing of a discrete-time (DT) improper-complex second-order cyclosta- 

^^ ' tionary (SOCS) random process is considered. In particular, it is of interest to find a pre-processing oper- 

ation that enables computationally efficient near-optimal post-processing. An invertible linear-conjugate 
linear (LCL) operator named the DT FREquency Shift (FRESH) properizer is first proposed. It is 
O ' shown that the DT FRESH properizer converts a DT improper-complex SOCS random process input 

to an equivalent DT proper-complex SOCS random process output by utilizing the information only 
^ , about the cycle period of the input. An invertible LCL block processing operator named the asymptotic 

r — ■ FRESH properizer is then proposed that mimics the operation of the DT FRESH properizer but processes 

L| ' a finite number of consecutive samples of a DT improper-complex SOCS random process. It is shown 

■«H- \ that the output of the asymptotic FRESH properizer is not proper but asymptotically proper and that its 

frequency-domain covariance matrix converges to a highly- structured block matrix with diagonal blocks 
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as the block size tends to infinity. Two representative estimation and detection problems are presented 
to demonstrate that asymptotically optimal low-complexity post-processors can be easily designed by 
exploiting these asymptotic second-order properties of the output of the asymptotic FRESH properizer 
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I. Introduction 

It is well known that the complex envelope of a real-valued bandpass wide-sense stationary 
(WSS) random process is proper, i.e., its complementary auto-covariance (a.k.a. the pseudo- 
covariance) function vanishes [H]. However, there are a lot of important improper-complex 
random processes that do not have vanishing complementary auto-covariance functions flU- 
H. For example, the complex envelope of a real-valued bandpass nonstationary signal is not 
necessarily proper and even the complex envelope of a real-valued bandpass WSS signal becomes 
improper in the presence of the imbalance between its in-phase and quadrature components [21. 
In order to fully capture the statistical properties of a complex-valued signal, including all the 
second-order statistics, the filtering of augmented signals has been proposed. This so-called 
widely linear (WL) filtering processes either the signal augmented by its complex conjugate or 
the real part of the signal augmented by the imaginary part [21, [111, where the former is referred 
to as the linear-conjugate linear (LCL) filtering [|7]. 

On the other hand, many digitally modulated signals are well modeled by wide-sense cyclo- 
stationary (WSCS) random processes [|8], [|9l, whose complex envelopes possess periodicity in 
all variables of the mean and the auto-covariance functions. To efficiently extract the correlation 
structure of a WSCS random process in the time and the frequency domains, the translation series 
representation (TSR) and the harmonic series representation (HSR) are proposed [[TOl . which are 
linear periodically time- varying (PTV) processings. Using these representations, it is shown [10] 
that a proper-complex WSCS scalar random process can be converted to an equivalent proper- 
complex WSS vector random process. Such second-order structure of a proper-complex WSCS 
random process has long been exploited in the design of many communications and signal 
processing systems including presence detectors [UTI . [[T2l . estimators [[T3l - [[T5l . and optimal 
transceivers [[T6l - [[T8l under various criteria. 

The complex envelopes of the majority of digitally modulated signals are proper and WSCS. 
However, as well documented in [[51, there still remain many other digitally modulated signals 
such as pulse amplitude modulation (R^.M), offset quaternary phase-shift keying (OQPSK), and 
Gaussian minimum shift keying (OMSK), that are not only improper-complex WSCS but also 
second-order cyclostationary (SOCS) [IT9l , i.e., the complementary auto-covariance function is 
also periodic in all its variables with the same period as the mean and the auto-covariance 
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functions. To exploit both cyclostationarity and impropriety, the LCL FREquency Shift (FRESH) 
filter has been proposed in [[20] that combines signal augmentation and linear PTV processing. 

Recently, another LCL PTV operator called the properizing FRESH (p-FRESH) vectorizer is 
proposed in [[T9l , [[211 by non-trivially extending the HSR. The p-FRESH vectorizer converts 
an improper-complex SOCS scalar random process to an equivalent proper-complex WSS vector 
random process by exploiting the frequency-domain correlation and complementary correlation 
structures that are rigorously examined in [Q. By successfully deriving the capacity of an SOCS 
Gaussian noise channel, it is shown that the optimal channel input is improper-complex SOCS in 
general when the interfering signal is improper-complex SOCS. It is also well demonstrated that 
such properization provides the advantage of enabling the adoption of the conventional signal 
processing techniques and algorithms that utilize only the correlation but not the complemen- 
tary correlation structure. These results warrant further research in communications and signal 
processing on the efficient construction and processing of improper-complex SOCS random 
processes. 

In this paper, we consider the block processing of a complex- valued random vector that is 
obtained by taking a finite number of consecutive samples from a discrete-time (DT) improper- 
complex SOCS random process. To proceed, the p-FRESH vectorization developed for continuous- 
time (CT) SOCS random processes is first extended to DT improper-complex SOCS random 
processes. Instead of straightforwardly modifying the CT p-FRESH vectorizer that properizes 
as well as vectorizes the input CT improper-complex SOCS random process, an invertible LCL 
operator is proposed in this paper that does not vectorize but only properizes the input DT 
improper-complex SOCS random process. Thus, the LCL operator is named the DT FRESH 
properizer and its output does not need vector processing. Specifically, the DT FRESH properizer 
transforms a DT improper-complex SOCS scalar random process into an equivalent proper- 
complex SOCS scalar random process with the cycle period that is twice the cycle period of the 
input process. 

To make this idea of pre-processing by properization better suited for digital signal processing, 
an invertible LCL block processing operator is then proposed that mimics the operation of the 
DT FRESH properizer. Although the augmentation of the improper-complex random vector 
by its complex conjugate can generate a sufficient statistic [[5l[, it is not only a redundant 
information of twice the length of the original observation vector but also improper. Although 
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the Strong uncorrelating transform (SUT) of the observation vector can generate a highly- 
structured sufficient statistic of the same length as the original observation vector [|21, [|22l . 
the output is still improper and, moreover, it requires the information about both the correlation 
and the complementary correlation matrices of the improper-complex random vector input for 
the transformation. Thus, neither the augmenting pre-processor nor the SUT allows the direct 
application of the conventional techniques and algorithms dedicated to the block processing of 
proper-complex random vectors. Motivated by how the DT FRESH properizer works in the 
frequency domain, the pre-processor proposed in this paper utilizes the centered discrete Fourier 
transform (DFT) and the information only about the cycle period of the DT improper-complex 
SOCS random process in order to convert a finite number of consecutive samples of the random 
process to an equivalent random vector. Unlike the output of the DT FRESH properizer, the output 
of this LCL block processing operator is not proper but approximately proper for sufficiently 
large block size. Thus, the pre-processor is named the asymptotic FRESH properizer. Specifically, 
the asymptotic FRESH properizer makes the sequence of the complementary covariance matrices 
of the output asymptotically equivalent f23\ to the sequence of all-zero matrices. 

In [[24|. it is shown that the complex- valued random vector consisting of a finite number of 
consecutive samples of a DT proper-complex SOCS random process has its frequency-domain 
covariance matrix that approaches a block matrix with diagonal blocks as the number of samples 
increases. This is because the periodicity in the second-order statistics of the DT SOCS random 
process naturally leads to a sequence of block Toeplitz covariance matrices as the number of 
samples increases, the sequence of block Toeplitz matrices is asymptotically equivalent to a 
sequence of block circulant matrices, and a block circulant matrix becomes a block matrix with 
diagonal blocks when pre- and post-multiplied by DFT and inverse DFT matrices, respectively. 
Since the asymptotic FRESH properizer mimics the DT FRESH properizer and the DT FRESH 
properizer outputs a DT proper-complex SOCS random process, it naturally becomes of interest 
to examine the asymptotic property of the frequency-domain covariance matrix of the asymptotic 
FRESH properizer output. It turns out that the output of the asymptotic FRESH properizer has 
the same property discovered in [|24]. 

Such properties of the covariance and the complementary covariance matrices may be used 
in designing, under various optimality criteria, low-complexity post-processors that follow the 
asymptotic FRESH properizer. In particular, a post-processor can be developed that approximates 
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the output of the asymptotic FRESH properizer by a proper-complex random vector having the 
block matrix with diagonal blocks that is asymptotically equivalent to the exact frequency-domain 
covariance matrix as its frequency-domain covariance matrix. Of course, this technique makes the 
post-processor suboptimal that performs most of the main operations in the frequency domain. As 
the two asymptotic properties strongly suggest, however, if the block size is large enough, then 
the performance degradation may be negligible. It is already shown in f24l that this is the case 
for the asymptotic property only of the covariance matrix, where a suboptimal frequency-domain 
equalizer is proposed that approximates the frequency-domain covariance matrix of a proper- 
complex SOCS interference by an asymptotically equivalent block matrix with diagonal blocks. 
It turns out that this equalizer not only achieves significantly lower computational complexity 
than the exact linear minimum mean-square error (LMMSE) equalizer by exploiting the block 
structure of the frequency-domain covariance matrix, but also is asymptotically optimal in the 
sense that its average mean-squared error (MSB) approaches that of the LMMSE equalizer as 
the number of samples tends to infinity. 

To demonstrate the simultaneous achievability of asymptotic optimality and low complexity by 
employing the post-processor that processes the output of the asymptotic FRESH properizer and 
exploits the two asymptotic properties, we consider two representative estimation and detection 
problems. First, for a DT improper-complex SOCS random signal in additive proper-complex 
white noise, a low-complexity signal estimator is proposed that is a linear function of the 
output of the asymptotic FRESH properizer. It is shown that the average MSE performance 
of the estimator approaches that of the widely linear minimum mean-squared error (WLMMSE) 
estimator as the number of samples tends to infinity. Second, for a DT improper-complex SOCS 
Gaussian random signal in additive proper-complex white Gaussian noise, a low-complexity 
signal presence detector is proposed whose test statistic is a quadratic function of the output of the 
asymptotic FRESH properizer It is shown that the test statistic converges to the exact likelihood 
ratio test (LRT) statistic that is a quadratic function of the augmented observation vector with 
probability one (w.p. 1) as the number of samples tends to infinity. Note that in both cases the 
asymptotic FRESH properizer as the pre-processor utilizes only the information about the cycle 
period. Thus, the adoption of adaptive estimation and detection algorithms may be possible that 
are developed for the processing of proper-complex random vectors. This advantage of employing 
the asymptotic FRESH properizer may be taken in other communications and signal processing 
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problems involving the block processing of DT improper-complex SOCS random processes. 

The rest of this paper is organized as follows. In Section II, definitions and lemmas related to 
SOCS random processes are provided and the DT FRESH properizer is proposed. In Section III, 
the asymptotic FRESH properizer is proposed and the second-order properties of its output are 
analyzed. In Sections IV and V, the application of the asymptotic FRESH properizer is considered 
to exemplary estimation and detection problems, respectively. Finally, concluding remarks are 
offered in Section VI. 

Throughout this paper, the operator E{-} denotes the expectation, the operator ^{x\n]} = 
X]^_oo x[n]e^'^'^'^^'^ denotes the discrete-time Fourier transform (DTFT) of an absolutely summa- 
ble sequence {x\n\)n, and the function 5{-) denotes the Dirac delta function. The sets Z and N 
are the sets of all integers and of all positive integers, respectively. The operator x denotes the 
Cartesian product between two sets. The superscripts *, ^, and ^ denote the complex conjugation, 
the transpose, and the Hermitian transpose, respectively. The operator * denotes the convolution!^ 
The operators and ® denote the Hadamard product and the Kronecker product, respectively. 
The matrices 1^, In, On, and Om,n denote the A^-by-A^ all-one matrix, the A^-by-A^ identity 
matrix, the A^-by-A^ all-zero matrix, and the M-by-A^ all-zero matrix, respectively. The matrix 
Pn denotes the A^-by-A^ backward identity matrix whose (m, n)th entry is given by 1 for 
m + n = N + I, and otherwise. The operator [A]m,n and tr{A} denote the (m,?7,)th entry 
and the trace of a matrix A, respectively. To describe the computational complexity, we will use 
the big-0 notation 0{g{N)) defined as f{N) = 0{g{N)) if and only if there exist a positive 
constant M and a real number A^o such that |/(A^)| < M|5f(A^)|, VA^ > A^o- 

II. DT FRESH Properizer 

In this section, the notion of DT second-order cyclostationarity is introduced and an LCL 
PTV operator is proposed that converts a DT improper-complex SOCS random process into 
an equivalent DT proper-complex SOCS random process. Similar to the p-FRESH vectorizer 
proposed in fil9] . where an input CT SOCS random process is converted to an equivalent 
CT proper-complex WSS vector random process, this operator as a pre-processor enables the 
adoption of the conventional signal processing techniques and algorithms that utilize only the 

* There should be no confusion from the superscript * that denotes the complex conjugation. 
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correlation but not the complementary correlation structure of the signal. Note that, unlike the 
p-FRESH vectorizer, this DT operator does not vectorize but only properizes the input improper- 
complex SOCS random process. A CT version of this properizer can be found in [|25l . 

A. DT SOCS Random Processes 

In this subsection, definitions and lemmas related to the DT SOCS random processes are 
provided. 

Definition 1: Given a DT complex- valued random process X[n] with a finite power, i.e., 
E{|X[r2]p} < oo,Vr2, the mean, the auto-correlation, the auto-covariance, the complementary 
auto-correlation, and the complementary auto-covariance functions of X [n] are defined, respec- 
tively, as 

fix[n]^E{X[n]}, (la) 

rx[n,m] = E{X[n]X[m]*}, (lb) 

cx[n,m] ^ E{iX[n] - fxx[n]){X[m] - fixHY}, (Ic) 

fx[n,m]=E{X[n]X[m]},aiid (Id) 

cx[n,m] ^ E{iX[n] - fix[n])iX[m] - f^xH)}. (le) 

Throughout this paper, all DT complex-valued random processes are assumed to be of finite 
power, i.e., Irxi^,""^]] < E{|X[n]p} < oo,Vn,Vm. 

Definition 2: The two-dimensional (2-D) power spectral density (PSD) Rxif, /') and the 2-D 
complementary PSD Rx{f, f) of a DT complex- valued random process X[n] are defined as 

oo oo 

RxUJ')=Y. E ^^K^]e"^'"^^"^^'"'^ (2a) 

m=— oo n=— oo 

and 

oo oo 

RxU.n=Y. E ^~^K"^]e-^'"^^'^-^''"\ (2b) 

m=— oo n=— oo 

respectively, if they exist. 

Since the 2-D PSD and the 2-D complementary PSD are the DT double Fourier transforms 
of rx[n^m] and fxln.m], respectively, they are always periodic in both variables / and /' with 
the common period 1. The set of all DT complex- valued random processes can be partitioned 
into two subsets by using the following definition. 
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Definition 3: [HI Definition 2] A DT complex- valued random process X[n] is proper if its 
complementary auto-covariance function vanishes, i.e., cx[n,m] = 0,V?7,,Vm, and is improper 
otherwise. 

Two types of stationarity can be defined as follows by using the second-order moments of a 
DT complex-valued random process. 

Definition 4: ll26l Section II-B] A DT complex- valued random process X[n] is second-order 
stationary (SOS) if, Vm, Vn, 

fJ'x[n] = Aix[0], (3a) 

cx[n,m] = cx[n — m,0], and (3b) 

Cx[n,m] = cx[n — 171,0]. (3c) 

Definition 5: A DT complex- valued random process X[n\ is SOCS with cycle period M E N 

if, Vn, Vm, 

l^x[n] = fix[n + M], (4a) 

cx [n,m] = cx[n + M,m + M], and (4b) 

cx[n, m] = cx[n + M,m + M]. (4c) 

We are mainly interested in the DT SOCS random processes in this paper. Note that a DT 
SOS random process can be viewed as a special DT SOCS random process with cycle period 1. 
Note also that the above definition of a DT SOCS random process is a straightforward extension 
of the definition of a CT SOCS random process in [[T9l . For the ease of comparison with the 
results in liT9l , we use the time indexes m and n in the orders appearing in Definitions [I]|5] 

In the following lemmas, the implications of the second-order cyclostationarity are provided 
in the time and the frequency domains, respectively. 

Lemma 1: For a DT SOCS random process X[n] with cycle period M E N, there exist 
irP[n])tL~o' and (rt^M)*!^ ^^ch that 

M~l 

rx [n, m] = Y, ^f [^ " ^je^ '"'="/''', (5a) 

fc=0 

and 

M-l 

rx [n, m] = J] ff \n - m](^^^'''^''"' . (5b) 

A:=0 
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Proof: Since r^ [n, I] = rx [n, n—l] is periodic in n with period M and is finite, there exist the 
DT Fourier series coefficients (?"x [^])fcio^ f^^ ^^'-h / such that r'-^[n,l] = ^^=^ rx[l]e^'^'^'^"'^^^- 
By replacing / with n — m, we obtain (|5al ). Similarly, we obtain (|5b] ). Therefore, the conclusion 
follows. □ 

Lemma 2: For a DT SOCS random process X[n] with cycle period Af E N, the 2-D PSD 
and the 2-D complementary PSD are given, respectively, by 



and 



-oo k=0 



oo A/-1 



(6a) 






i=— oo fc=0 



(6b) 



where R^if) = ^{rP[n]} and ^5J^^(/) = ^{rP[n]}. 

Proof: It is straightforward by applying the definitions (|2a1 ) and (|2bl) in Definition [2l 
respectively, to (l5al) and (|5bl) in Lemma [T] □ 

Note from (l6]) that both Rxif, f) and Rx{f, /') consist of 1/M-spaced impulse fences along 
the lines f = f ± k/M, \fk E Z. The above lemmas are the extensions of the results in [[T9l for 
a CT SOCS random process. 

B. DT FRESH Properizer 

In this subsection, an LCL PTV operator is proposed that always outputs an equivalent DT 
proper-complex random process, regardless of the propriety of the input DT SOCS random pro- 
cess. In what follows, all SOCS random processes are DT processes unless otherwise specified. 

To proceed, the following frequency-selective filter is defined that has half of the entire 
frequency band as its stopband. 

Definition 6: Given a reference frequency 1/M, a linear time-invariant system with impulse 
response gM\n\ is called the frequency-domain raised square wave (FD-RSW) filter if its fre- 
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Fig. 1. Frequency response of the FD-RSW filter with reference frequency 1/M. 



quency response G'a/(/) = •^{gM[n]} is given by 



GmU) 



0, 


for -^</<0, 


1, 


for < / < 2^7, and 


Gm if + jj), 


elsewhere. 



(V) 



Fig. [H shows the frequency response of the FD-RSW fiher with reference frequency 1/M, 
which akemately passes the frequency components of the input signal in every other interval 
of bandwidth 1/(2A/). Hereafter, Qj\j denotes the support of this FD-RSW filter. By using the 
impulse response QMln] of the FD-RSW filter with reference frequency 1/M, we can define an 
LCL PTV filter as follows. 

Definition 7: Given a reference frequency 1/Af and an input X[n], the DT FRESH properizer 
is defined as a single-input single-output LCL PTV system, whose output is given by 

Y[n] ^ X[n] * toN + {X[n]* * gM[n]) e^^WCsAf). (g) 

Fig. [2l shows how the DT FRESH properizer works in the time domain, where X[n] is the 
input, Y[n] is the output, and Xi[n] and X2[n\ are, respectively, the first and the second terms 
on the right side of dS). Note that, unlike the upper branch where X[n] is processed by the 
FD-RSW filter to generate Xi[r2], X[n]* is processed in the lower branch by the FD-RSW filter 
and multiplied by a complex-exponential function e^J^'^"/'^^*'^^ to generate X2[n]. Even though 
— l/(2Af) is chosen as the frequency of the complex-exponential function, any (2A; + 1)/(2M) for 
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Fig. 2. DT FRESH properizer with reference frequency 1/AI viewed in tfie time domain. 



- Y[n] 



k E Z can be chosen. The reason why this is so becomes clear once the DT FRESH properization 
is viewed in the frequency domain. 

Fig. [3] shows how the DT FRESH properizer works in the frequency domain, especially when 
the input is a deterministic signal s[n] with the DTFT S{f) = ^{s[n]}, the outputs of the upper 
and the lower branches in Fig. [2] are denoted by si[n] with the DTFT 5'i(/) = =^{si[n]} and 
S2[n\ with the DTFT 5*2 (/) = ^{s2[n]}, respectively, and the output signal is denoted by t[n] 
with the DTFT T(/) = ^{t[n]}. Note that S{f) is processed by the FD-RSW filter to generate 
the first term Si{f) of the DTFT of the DT FRESH properizer output, while S{—f)* is processed 
by the FD-RSW filter and shifted in the frequency domain to generate the second term 5'2(/). 
Thus, 5*1 (/) contains all the frequency components of s[n] on the support ^ of the FD-RSW 
filter, while 5*2 (/) contains all the remaining frequency components. Since the supports of 5'i(/) 
and 5*2 (/) do not overlap, the DTFT of the output T(/) of the DT FRESH properizer contains 
all the frequency components of the input signal S{f) without any distortion. Note also that, 
due to the periodicity of the DTFT with period 1, the frequency shift of 5*2 (/) by any k/M 
for /c G Z generates the output signal that contains the same information as the input does. The 
following lemma makes this invertibility argument more precise. 

Lemma 3: From the output Y[n] of the DT FRESH properizer with reference frequency 1/M, 
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Fig. 3. Example that shows how the DT FRESH properizer with reference frequency 1/M works in the frequency domain, 
when the input is a deterministic signal s[n\. 



the input X[n\ of the DT FRESH properizer can be recovered as 



X[n] = Y[n]*gM[n] + {{Y[n]e^'--/^'''^) * gMln]}' 



(9) 



Proof: Straightforward by substituting (|8]) into the right side of (I9l). □ 

Note that, at this point, the DT FRESH properizer may not be more than one of many possible 

invertible operators. The reason why this operator is named the properizer will become clear 
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once the second-order property of the output is analyzed as follows when its input is a zero-mean 
SOCS random process. 

Theorem 1: If the input X[n] to the DT FRESH properizer with reference frequency 1/M 
is a zero-mean SOCS random process with cycle period M, then the output Y[n] becomes a 
zero-mean proper-complex SOCS random process with cycle period 2M, i.e., the mean, the 
auto-correlation, and the complementary auto-correlation functions of Y[n] satisfy 

/iy[n] = E{r[n]} = 0, (10a) 

rY[n,m] = E{Y[n]Y[m]*} = rY[n+2M,m + 2M], and (10b) 

fY[n,m] = E{Y[n]Y[m]} = 0, (10c) 

Vm,Vr2, respectively. 

Proof: It is straightforward to show (llOal ) by using /ix[^] = 0,Vri. Let Xi[n] and X2[n] be 
defined again as shown in Fig. [21 Then, the auto-correlation function ry[?7,,m] of Y[n] can be 
written as 

rY[n,m] = E{Xi[n]Xi[m]*} + E{Xi[n]X2[m]*} + E{X2[n]Xi[m]*} + E{X2[n]X2[m]*}. (11) 

The first term on the right side of (fTTT i can be rewritten as 

E{X,[n]X,[m\*} = [ [ GM{f)e^^^f''Rx{fJ')GM{rre-^''^^''^dfdf' (12a) 

JO Jo 

= E ( / Rx if - Y^ ^'^'^"^df) e^^-*-, (12b) 

fc=0 ^-^'^ ^ ^ ^ 

where (I12al) holds by Parseval's relation and (1 12b I ) holds by substituting (l6al) into (I12al) . It turns 

out in (I12bl) that E{Xi[r2]Xi[m]*} is periodic in both n and m with period M. Similarly, the 

second term of (fTT)) can be rewritten as 

E{XMX2[mr} = E ( / ^x (/ - ^) e^^-^^— )rf/) e^^^^-. (13) 

It also turns out in (fT3l) that E{Xi[?2]X2[m]*} is periodic in both n and m with period 2M. 
In the same way, the other two terms can be obtained, which turn out to be periodic in n and 
m with period 2M and M, respectively. Thus, the auto-correlation function rY\n,m\ of Y[n] 
satisfies (llObI) . 
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On the other hand, the complementary auto-correlation function fY[n,m] of Y[n] can be 
written as 

fY[n,m] = E{Xi[n]Xi[m]} + E{Xi[n]X2[m]} + E{X2[n]Xi[m]} + E{X2[n]X2[m]}} . (14) 

The first and the second terms of the right side of (fT4] ) can be rewritten, respectively, as 

E{X,[n]X,[m]}=[[ RxU,-fy^^^''^'^^''''dfdf (15a) 



and 



E{X^[n]X2[m]}=ff Rx{f, -neMf-+f'-^-^)dfdf', (15b) 

7 J^xSf 



Vn, Vm by using Parseval's relation. These two terms are all zeros because the impulse fences of 
Rxif, —f) and Rxif, —f) along the line / = —f'~k/M for any k do not cross the integration 
area ^ x ^. Fig. |4] illustrates these lines and the integration area. Similarly, the other two terms 
of fY[n,m] can be obtained, which again turn out to be all zeros. Thus, the complementary 
auto-correlation function fY[n,m\ of Y[n] satisfies (llOcI ). Therefore, the conclusion follows. □ 

This theorem shows that the DT FRESH properizer in general doubles the cycle period at the 
cost of the propriety of the output. However, it does not double the cycle period if the input is 
already proper. 

Corollary 1: If the input X[n\ to the DT FRESH properizer with reference frequency 1/M is 
a zero-mean proper-complex SOCS random process with cycle period M, then the output Y[n] 
is a zero-mean proper-complex SOCS random process with cycle period M. 

Proof: It suffices to show that ry [n, m] is periodic in n and m with period M. It is already 
shown in Theorem [H that the first and the forth terms of rY[n, m] are periodic in n and m with 
period M. As it can be easily seen in (fT3l) . the second term E{Xi[?7,]X2[m]*} on the right side 
of ([TT]) is zero, Vn, Vm, because the propriety of X[n] implies Rxif, f) = 0, V/, V/'. Similarly, 
the third term is zero. Therefore, the conclusion follows. □ 

It is already shown that the amount of the frequency shift in the second term of ^ can be 
any (2k + 1)/(2M), for k E Z,to satisfy the invertibility of the DT FRESH properizer. It can be 
also shown that Theorem [T] holds for any frequency shift {2k + 1)/(2M), for A; G Z. Moreover, 
since any integer multiple of M is also a cycle period of an SOCS random process with cycle 
period M, the random process can be FRESH properized by using any reference frequency 
l/{kM),yk E N. Thus, it is not unique to FRESH properize an SOCS random process. 
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Fig. 4. Solid lines represent the impulse fences of Rx{f,—f'), Rx{—f,f'), Rx(f,—f'), or Rx{—f,f')- Shaded area 
represents the integration area ^ x ?#. 



III. Asymptotic FRESH Properizer 



In this section, the block processing of an SOCS random process is considered. Motivated 
by how the DT FRESH properizer works in the frequency domain, an LCL block operator is 
proposed that converts a finite number of consecutive samples of an SOCS random process to an 
equivalent random vector. Unlike the DT FRESH properizer proposed in the previous section, 
this invertible operator does not directly make the complementary co variance matrix of the output 
vector vanish. Instead, it is shown that the LCL operator makes the complementary covariance 
matrix of the output vector approach all-zero matrix as the number of samples tends to infinity. 
This is why it is named the asymptotic FRESH properizer. 
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A. Asymptotic FRESH Properizer and Its Inverse Operator 

Let X be the length-MA^ vector obtained by taking the MN consecutive samples of a DT 
signal, where and in what follows it is assumed that A^ is a positive even number. Motivated 
by the DT FRESH properizer, we introduce in this subsection an LCL block operator and its 
inverse. 

To proceed, some definitions are provided. 

Definition 8: The centered DFT matrix Wmn is defined as an {MN)-hy-{MN) matrix whose 
(m, ?i)th entry, for m, n G {1, ■ ■ ■ , MN}, is given by 

A 1 



\W, 



MN\m,n 



-j27r(m-CMJv)(n-Cj\/iv)/(A-fA^) 



(16) 



Vmn 

where cmn = (MN + l)/2. 

It is well known that the matrix-vector multiplication with a centered DFT matrix can be 
implemented with low computational complexity [|27l , as the multiplication with an ordinary 
DFT matrix is efficiently implemented by using the fast Fourier transform algorithm. 

Definition 9: Given M and an even number N, the {MN)-by-(MN) matrix Gm,n is defined 
as 



G 



M,N 



IM 



On/2 On/2 
On/2 In/2 



(17) 



Similar to the FD-RSW pulse, the matrix Gm,n will be called the raised square wave (RSW) 
matrix because, when pre-multiplied to a column vector or a matrix, it turns the {{m — l)N + n)th 
row, for m = 1, ■ ■ ■ , M and n = 1, ■ ■ ■ , N/2, into all zeros, i.e., it alternately nulls every other 
band of N/2 consecutive rows. 

Definition 10: Given M and an even number N, the {MN)-by-{MN) matrix Sm,n is defined 
as 



>M,Af 



o 



MN-N/2,N/2 



'-N/2 



'■MN-N/2 



o 



N/2,MN-N/2 



(18) 



Note that the matrix Sm.n, when pre-multiplied to a column vector or a matrix, circularly 
shifts the rows by N/2, which corresponds to multiplying e^J^'^"/'^^^^) in the second term of the 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY (VER 1.0) 17 

right side of ([8]). Now, we are ready to introduce an LCL operator that is the block-processing 
counterpart to the DT FRESH propertizer. 

Definition 11: Given M and an even number N, the LCL operator / with input x and output 
y = f{x), both of length MN, is called the asymptotic FRESH properizer if 

fix) ^ W]Ji^{Gm,nWmnx + Sm,nGm,nWmnX*). (19) 

Note that the input x is pre-multiplied by the centered DFT matrix Wmn and the RSW matrix 
Gm,n, while the complex conjugate x* is multiplied additionally by the circular shift matrix 
Sm,n to generate the frequency-domain output WmnV- 

Fig. \5\ shows Wmnx, Gm,nWmnx, and Sm,n Gm,nWmnX*, when the /th entry of the 
WmnX is denoted by x^. Note that, similar to the DT FRESH properization illustrated in Fig. [3l 
the locations of all possible non-zero rows of Gm,nWmnx and Sm,nGm,nWmnx* do not 
overlap, which makes WmnV contain all the entries of Wmnx without any distortion. Note 
also that the amount of the circular shift of Sm,nGm,nWmnX* by any kN for k E Z generates 
the signal that contains the same information as Sm,nGm,nWmnx* does. The following lemma 
makes this invertibility argument more precise. 

Lemma 4: From the output y of the asymptotic FRESH properizer with parameters M and 
A^, the input x of the asymptotic FRESH properizer can be recovered as 

X = Wl^jj^{GM,NWMNy + PMNGM^NSlj^^Wljj^y*) (20a) 

= r\y) (20b) 

Proof: By substituting (fT9l) into the right side of (I20al) . we have WIjj^GIjj^Wmnx + 
^MN^M,N Sm,nGm,nWmnx* + W^j^PmnGm,nS^j^^Gm,n W&in^* + W^jj^PmnGm,n- 
Sjj j^Sm,nGm,nW^jj^x. The second and the third terms vanish since Gm,nSm,nGm,n = Omn 
and Gm,nS1i n^m,n = Omn, respectively. Thus, the right side of (I20al l becomes W2/jv(G^^ jv + 
PmnGm^j^Pmn)Wmnx. Moreover, it can be easily shown that G|^jv = Gm,n and Gm,n 
+PmnGm,nPmn = Imn, because Pmn = Pm ® Pn by the properties of the Kronecker 
product. Therefore, the conclusion follows. □ 

Note that, at this point, the asymptotic FRESH properizer may not be more than one of many 
possible invertible operators. The reason why this operator is named the asymptotic properizer 
will become clear once the second-order property of the output is analyzed in the next subsection 
when its input is a finite number of consecutive samples of a zero-mean SOCS random process. 
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Fig. 5. Illustration that shows how the asymptotic FRESH properizer works in the frequency domain. 



B. Second-Order Properties of Output of Asymptotic FRESH Properizer 



Let X be the length-2Af A^ augmented vector defined as 



- A 
X = 



X 



x^ 



(21) 
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where and in what follows it is assumed that x consists of a finite number of consecutive 
samples of a zero-mean SOCS random process X[n\ with cycle period M G N. Then, the output 
y = f{x) in (fT9l ) of the asymptotic FRESH properizer can be rewritten as 



y = WI^j^Gm,nWmnx 



(22a) 
(22b) 



where the (M7V)-by-(2MiV) matrix Gm,n and the (2MA^)-by-(2MA^) matrix Wmn are given 
by 



G 



M,N 



Wmn = 



Gm,n Sm,nGm,n 

Wmn Omn 
Omn Wmn 



and 



(23a) 



(23b) 



respectively, and y is the centered DFT of y. Thus, the covariance matrix Ry = E{i/y^} and 
the complementary covariance matrix Ry = E{j/?/^} of y are given by 



-R,, 



Wmn^m,nWmnRxWmn^m nWmn 



and 



rV. 



tT niT 



Ry — WmnGm,nWmnRxWmnGm nWIiis/, 



(24a) 



(24b) 



respectively, where R^ = E{a;x^} and Rg. = 'E{xx^} denote the covariance and the comple- 
mentary covariance matrices of the augmented vector x, respectively. 
Throughout this paper, we call 



Ry ^ E{yy^} = WmnRvWI^ 



MN 



and 



Ry^E{yy'} = WMNRyW, 



T 

MN 



(25a) 



(25b) 



the frequency-domain covariance and the frequency-domain complementary covariance matrices, 
respectively. 

To analyze the asymptotic second-order properties of y, we briefly review definitions and 
related lemmas for asymptotic equivalence between two sequences of matrices. 
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Definition 12: [|23l Let {Ak)k and (3^)^ be two sequences of N^-hy-Nk matrices with 
A'^fc, — 7- oo as A; — 7- oo. Then, {Ak)k and {Bk)k are asymptotically equivalent and denoted by 
Ak ~ -Bfc if 1) the strong norms of A^ and Bf^ are uniformly bounded, i.e., there exists a constant 
c such that || A^H, ||-Bfc|| < c < oo, VA;, and 2) the weak norm of A^ — Bk vanishes asymptotically, 
i.e., linifc^oo \Ak — Bk\ = 0, where the strong norm ||A|| and the weak norm |A| of an A^^- 



by-A^fc matrix A are defined as ||A|| = maxa^^o y x^ A^ Ax / x^x and |A| = -^ti{AA'^)/Nk, 
respectively. 

Lemma 5: If two sequences {Ak)k and {Bk)k oi Nk-hy-Nk matrices are asymptotically equiv- 
alent and if the strong norms of N^-hy-Nk matrices Ck and D^ are uniformly bounded, then 

CkAkDk ~ CkBkDk. 

Proof: Since Ck ~ Ck and Dk ~ Dk, it is straightforward to show the conclusion by 
applying |[23l Theorem l-(3)] twice. □ 



Lemma 6: If Ak ~ Bk and Ck ~ Dk, then (A^ + Cfc) ~ (-6,^. + Dk). 

Proof: By applying the triangle inequality, it is straightforward to show that the strong norms 
of Ak + Ck and Bk + Dk are uniformly bounded and that the weak norm of Ak + C^ — Bk — Dk 
vanishes asymptotically. Therefore, the conclusion follows. □ 

For brevity, in what follows, two square matrices Ak and Bk of the same size will be said 
to be asymptotically equivalent if Ak ~ Bk- Now, we introduce our definition of asymptotic 
propriety. 

Definition 13: A sequence {yk)k of length-Nk complex-valued random vectors is asymptoti- 
cally proper if the complementary covariance matrix of yk is asymptotically equivalent to the 
Nk-by-Nk all-zero matrix. 

For brevity, a random vector yk will be said to be asymptotically proper if the sequence {yk)k 
is asymptotically proper. Before showing that the output y of the asymptotic FRESH properizer 
is asymptotically proper, the second-order properties of the augmented vector x are examined 
as follows. 

Proposition 1: If the random vector x is obtained by taking the MN consecutive samples 
from a zero-mean SOCS random process with cycle period M, then the covariance matrix Rx 
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and the complementary covariance matrix R^ of the augmented vector x defined in (|2TI) satisfy 



WmnRo^WI^n ~ " (26a) 



and 



WmnR^WZn ~ ^, 



respectively, where the (2MA^)-by-(2MA^) matrices ft and Ct are given by 

n ^ (WmnR^WJ^jn) (l2M ® In) 



and 



O ^ (V^MiV^^VFA';^) (l2M fV), 



(26b) 



(27a) 



(27b) 



respectively. 

Proof: By definition, the left sides of (I26al) and (I26bl) can be rewritten, respectively, as 



WmnR^W^in 



WmnR^W^in WmnR^W^in 



WmnRIW^ 



MN 



WmnRIW^ 



MN 



(28a) 



and 



WmnR^WIjj, 



WMMR.WTr^r WmnR^WZ 



MN-Ttxyv MN 



MN-I^x '''^ MN 



wmnr:wi,j, wmnr:wi,j. 



(28b) 



Since both R^ and R^ are {MN)-hy-{MN) block Toeplitz matrices with block size M-by-M, 
each submatrix on the right side of (12 8 al l is asymptotically equivalent to the Hadamard product 
of {1m ® In) and the submatrix itself as shown in ll24l Proposition 1]. Thus, by using the fact 
that A ~ B if each submatrix of A is asymptotically equivalent to the corresponding submatrix 
of B [ [24l Lemma 5], we obtain (|26al) . Similar to (I26al) . it can be shown by using the property 
Wmn = ^mnRmn of the centered DFT matrix that each submatrix on the right side of (I28bl) 
is asymptotically equivalent to Hadamard product of (1m ® Pn) and the submatrix itself. Thus, 
we obtain (I26bl) . Therefore, the conclusion follows. □ 
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Fig. 6. Magnitudes of exemplary (a) WmnR^WIJj^ for iV = 50, (b) 17 for iV = 50, (c) WmnR^Wj^jj^ for iV = 500, 
and (d) n for iV = 500, when M = 2. 



Figs. [6] and |7] show the entry-by-entry magnitudes of the exemplary frequency-domain co- 
variance and complementary covariance matrices. To generate a zero-mean improper-complex 
SOCS random process, equally-likely BPSK symbols are linearly modulated with the square- 
root raised cosine (SRRC) pulse having roll-off factor 0.5 and transmitted at the symbol rate 
of I/T5 [symbols/sec]. After the BPSK signal passes through the frequency selective channel 
with impulse response h{t) = 0.22e^''/^6{t - 5T,) + 0.44eJ^/2^(t - 3T,) + 0.87eJ^/*^5(t - T,), 
the SOCS random process X[n\ is obtained by 2-times over-sampling the received signal, i.e., 
the sampling rate is 2/Ts [samples/sec]. Thus, the zero-mean SOCS random process X[n] has 
the cycle period 2. As shown in (|26al) . it can be seen in Fig. [6] that WmnRxW^j^ approaches 
CI as N increases. Similarly, as shown in (I26bl l. it can be seen in Fig. |7] that WmnRxW^^j^ 
approaches f2 as A^ increases. 
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Fig. 7. Magnitudes of exemplary (a) WMNRa^Wjit^ for N = 50, (b) O for iV = 50, (c) V^Miv-RsW^A/iv for N = 500, 
and (d) fi for iV = 500, when U = 2. 



Now, the asymptotic second-order properties of the output y of the asymptotic FRESH 
properizer are examined as follows. 

Theorem 2: If the input x to the asymptotic FRESH properizer with parameters M and A^ 
is obtained by taking the MN consecutive samples from an SOCS random process with cycle 
period M, then the covariance and the complementary covariance matrices of the output y of 
the asymptotic FRESH properizer satisfy 



Ry ~ W'^.^T.Wmn 



(29a) 



and 



Ry ~ OmN) 



(29b) 
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where an {M N)-hy-{M N) matrix S is given by 



S — Gm,n^G]^^j^ (30a) 

= RyQ{l2M(S)lN/2). (30b) 

Proof: By applying Lemmas [5] and [6] to (I24al) and (I26al) . we obtain i?y ~ Gm,n^G'^i n- 
Since Gm,n^GPk4,n can be rewritten as (GA/,Arii^MAr-RsW^l/ArG^f,jv)0(GAf,Ar(l2M®-f7v)<5^/,7v) 
by the definition (|27al) . we obtain ii^^ ~ S from Gm,n{'^2M®In)G'^,^ j^=12m®In/2- Similarly, 
Ry ~ (Ga/,7vW^a/^^s1^a^^GI,,v) © (Gm,7v(12M ® Pn)GIj^j,) by dalS), (I26bl), and dlTbl). 
Thus, we obtain Ry ~ Omn from GA/,Ar(l2M ® Pn)G1,i n = ^a/tv- Therefore, the conclusion 
follows by Lemma [51 □ 

By Definition [13] and the above theorem, the output y of the asymptotic FRESH properizer is 
indeed asymptotically proper. Fig. [8]-(a) illustrates how the frequency-domain covariance matrix 
Ry of the output y of the asymptotic FRESH properizer is constructed from WmnRx^mn- 
The thick diagonal lines in each block of size A^-by-A^ represent all possible non-zero entries 
of 17 that is asymptotically equivalent to WmnRsWmn- Then, Ry that is asymptotically 
equivalent to S is obtained by collecting the shaded sub-blocks of size (A^/2)-by-(iV/2). Fig. [8]- 
(b) illustrates how the frequency-domain complementary covariance matrix Ry of the output 
y of the asymptotic FRESH properizer is constructed from WmnRx^mn- The thick anti- 
diagonal lines in each block of size A^-by-A^ represent all possible non-zero entries of f2 that 
is asymptotically equivalent to WmnRxWIi^^. Then, Ry that is asymptotically equivalent to 
Omn is obtained by collecting the shaded sub-blocks of size (A^/2)-by-(A^/2). 

It is already shown that the amount of the circular shift of Gm^n^mn^* in ([TVl) can be 
any N{2k + l)/2, for A; G Z, to satisfy the invertibility of the asymptotic FRESH properizer. It 
can be also shown that Theorem [2] holds for every circular shift by A^(2A; + l)/2, for k E Z. 
Moreover, since Rx and R^ that are (MA^)-by-(MA^) block Toeplitz matrices with block size 
M-by-M can be viewed as block Toeplitz matrices with block size (/M)-by-(/M), V/ G N, the 
sampled vector x can be asymptotically properized by using any parameters IM and N/l, for 
all / G N such that N/l is even. Thus, similar to the DT FRESH properizer, it is not unique to 
asymptotically properize a finite number of consecutive samples of a zero-mean SOCS random 
process with cycle period M. 
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Fig. 8. Illustrations that show how to construct Ry and Ry from WmnRs^WI'ji^ ~ fi and WmnRs^Wmn ~ ^' 
respectively. Thick diagonal and anti-diagonal lines represent all possible non-zero entries of $7 and f2, respectively. 



In the following two sections, two representative estimation and detection problems are pre- 
sented to demonstrate that low-complexity asymptotically optimal post-processors can be easily 
designed by exploiting the asymptotic second-order properties of the output of the asymptotic 
FRESH properizer. 



IV. Application of Asymptotic FRESH Properizer to Signal Estimation Problem 

In this section, given a finite number of consecutive samples of a zero-mean improper-complex 
SOCS random process in additive proper-complex white noise, the minimum mean squared error 
estimation of the random process is considered when only the second-order statistics of the 
observation vector are provided. It is well known [|5] that an optimal estimator is the WLMMSE 
estimator that linearly processes the observation vector augmented by its complex conjugate. This 
optimality of the use of the augmented vector still holds even if the output of the asymptotic 
FRESH properizer is processed, because the output as an equivalent observation vector is still 
improper. Instead, motivated by Theorem |2l we propose to use the LMMSE estimator that 
processes the output of the asymptotic FRESH properizer as if the signal component is a 
proper-complex random vector with its frequency-domain covariance matrix being equal to the 
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masked version of the exact frequency-domain covariance matrix. It turns out that this suboptimal 
linear estimator is asymptotically optimal in the sense that the difference of its average MSE 
performance from that of the WLMMSE estimator having higher computational complexity 
converges to zero as the number of samples tends to infinity. 

A. Asymptotically Optimal Low-Complexity Estimator 
Let r be the length-MA^ observation vector modeled by 

r = X + V, (31) 

where x is the desired signal to be estimated that consists of the MN consecutive samples of 
a zero-mean improper-complex SOCS random process with cycle period M E N, and v is the 
additive proper-complex white noise vector. 

Given the observation model (ISTT ) and the second-order moments E{x} = Omn, E{xx^} = 
Rx, 'E{xx'^} = Rx, E{v} = Omn, E{vv'^} = a'^I^N, and 'E{vv'^} = Omn of x and v, the 
WLMMSE estimation problem can be formulated as 

minimize e^^' 

a'wl (32) 

subject to x^i = Fir + F2r*, 
where the average MSE e^^^ is defined as 

e^''^^j^E{\\x-x^,r}. (33) 

In the following lemma, the optimal solution to the above estimation problem is provided along 
with its average MSE performance. It is well known [5] that the optimal estimator that minimizes 
the average MSE without the wide linearity constraint becomes this WLMMSE estimator if the 
desired signal vector and the additive noise vector are both Gaussian. 

Lemma 7: The WLMMSE estimator a^opt of x as the solution to the optimization problem 
(|321) is given by 

iCopt = i^i,opt^ + i^2,optr'* (34) 

where 

-F'l,opt = {Rx — RxRr R*x){Rr — RxRr Rx) (35a) 

and 
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2, opt 



yRx ~ Rx-^r -^xji-^r ~ -^x-^r ^^x) ■ (35b) 



with Rr = E{rr^} = a^lAjN+Rx- Moreover, the average MSE e^^ - E{\\x-Xoptf}/{MN) 
is given by 

^opt = "aJ/y ^^(-^^ ~ Fi,optRx - F2^optRx)- (36) 

Proof: See [51 Section 5.4] and the references therein. □ 

The WLMMSE estimator in (|34|) requires the computation of two matrices -Fi,opt and -F2,opt 
to be pre-multiplied to r and r*, respectively. As it can be seen from (13 5 al l and (I35bl ). the 
major burden in computing these matrices comes from the multiplications and the inversions of 
(M N) -by -(MN) matrices. It can be shown that the matrix multiplications in computing -Fi,opt 
and F2opt require 0{M^N^) complex-valued scalar multiplications because the matrix R~^ 
does not have any structure to be exploited in matrix-matrix multiplications even though Rr 
is block Toeplitz. By the same reason, the matrix inversion {Rr — RxRr^* R^) in (I35al) and 
(I35bl) requires 0{M^N^) complex-valued scalar multiplications. Thus, the overall computational 
complexity of the WLMMSE estimator is C(M^iV^). 

An altemative observation model that is equivalent to the original one in (ISTT ) can be obtained 
by applying the asymptotic FRESH properization to r as 

s = y + w, (37) 

where s = f{r), y = f{x), and w = f{v). Since the equivalent observation vector s is not 
proper in general, the linear processing of s and s* is still needed for the second-order optimal 
estimation of x. Instead, we propose to use a suboptimal linear estimator as follows, which is 
a function only of s. 

Definition 14: The proposed estimate Xp of x is defined by 

Xp = r\FpS), (38) 

where Fp is given by 

Fp ^ Wlij.'^ia^lMN + ^)-'Wmn. (39) 

Note from dM) and ^ that i:{a^lMN + ^y^WAfNS is the estimate of WmnU = WMNfix) 
and, consequently, that FpS is the estimate of y = f{x). Let yp = FpS. Then, the estimate Xp 
of X = f~^{y) is obtained as (l38l) by applying the inverse operation f ^ to yp. 
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It can be immediately seen that Fp is chosen to make yp the LMMSE estimate of y when the 
covariance and the complementary covariance matrices of y are set equal to W^j^YIWmn and 
Omn appearing in (|29a| ) and (|29bk respectively. This idea of using the asymptotically equivalent 
matrices is motivated by Theorem |2l which naturally leads to the asymptotic optimality of the 
proposed estimator as shown in what follows. 

First, the invariance of the Euclidean norm under the asymptotic FRESH properization is 
shown. 

Lemma 8: The input x and the output y = f{x) of the asymptotic FRESH properizer have 
the same Euclidean norm, i.e., ||x|| = \\y\\. 

Proof: Since the centered DFT matrix Wmn is unitary, we have ||a:;|| = || Wj\/7va3||. Let 
yi and yl denote the /th entries of y and Wmnx, respectively. Then, we have ||WA/Ara:;p = 
^1=1 lz//P = Z]«=i Ivil"^ = II2/P by Definition \TT\ of f{x) in (fT9l) . which is illustrated in Fig. |5] 
Therefore, the conclusion follows. □ 

Second, the average MSE of the proposed estimator is derived. 

Lemma 9: The average MSE Sp = E{||a3 — a:;pp}/(Af A^) of the proposed estimator is given 
by 

4""^ = ^ tr{S - i:{a'lMN + S)-iS}. (40) 

Proof: The linearity of the asymptotic FRESH properizer leads to y — yp = f{x — Xp). Then, 
by Lemma [8l the average MSE performance can be written as Ep = 'E{\\y — yp\\'^}/{MN). 
Since Wmn is unitary and yp = FpS, it can be shown that we have 

4"^^ = ]^tr{K^ - 2^ia'lMN + ^r'Ry 

^^{a^lMN + T.)-\a^lMN + Ry){o^lMN + T^Y^T.) . (41) 

It can be also shown that if A = A {I2M ® IN/2) and B = B © (I2A/ ® IN/2) for {MN)- 
hy-{MN) matrices A and B then tr(A) = tr(A) and tr(AB) = ti{AB). Thus, the right side 
of (|4TI) is invariant under replacing R^ with S = R^ © {I2M ® In 12) in (I30bl) . Therefore, (1411 
can be simplified to (|40|) . □ 

Now, the asymptotic optimality of the proposed estimator is provided. 

Theorem 3: The average MSE of the proposed estimator approaches that of the WLMMSE 
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estimator as the number of samples tends to infinity in the sense that 

i.m(4"'-C>)-o- («) 

Proof: It is shown in [|5l Section 5.4] that the average MSE e:^} in ^^ of the WLMMSE 
estimator can be simplified as 

£^opt = 9/v#/v ^^'l-^^ ~ ^^{^ I2MN + Rx) Rx}- (43) 

Now, we define Q = Wj^jj^CIWmn and introduce 

4S = ^ tr{f^ - nia'hMN + 0)-i^} (44a) 

= ^^ tT{n - Qia'hMN + f^)-'^}, (44b) 

which is obtained by replacing Rx in (l43l) with S7 and by using the fact that Wmn is unitary. 
Note that WmnRxW^J^ ~ fi as shown in (I26al) . This introduction of fopt is motivated by 
the result in [[24l Theorem 1], where a suboptimal estimator is proposed for the estimation of a 
proper-complex cyclostationary random signal and its average MSE is shown to approach that 
of the LMMSE estimator as the number of samples tends to infinity. Since Rx is a covariance 
matrix, it is positive semidefinite, which implies that all the eigenvalues of (cr^/2Miv + Rx)~^ 
in P3l) are upper bounded by l/a'^. Note that fl and tl are also positive semidefinite because 
(I2M ® IN/2) in (I27al) is positive semidefinite. Similarly, all the eigenvalues of {o-'^I2mn + fi)^^ 
in (I44al) are also upper bounded by l/a^. Thus, the strong norms ||(cr^/2j\/iv + Rx)^^]] and 
||(cr^/2AfAr + ri)~^|| are uniformly upper bounded by l/o-^ for any matrix size. Recall that if 
Ak ~ -Bfc and ||A^^||, ll-B^"*^!! < c < 00, V/c, for a positive constant c, then A^^ ~ Bj^^ 
[|23l Theorem l-(4)]. Thus, we have {a'^hMN + -R*)"^ ~ {o-'^hMN + ri)-^ Recall also that 
if two sequences {Ak)k and {B^jk of A^fc-by-A'^fc matrices are asymptotically equivalent then 
lim^^oo tr(Afc — Sfc)/A''fc = [|23l Corollary 1]. Thus, combined with Lemmas |5] and |6l we have 
limAr_i.oo(£^opt ~^opt ) = 0- Therefore, in order to show (|42T i. it now suffices to show i^p^ = Sp . 
Let E2M,N be the (2MA^)-by-(2MA^) matrix that permutes the rows of the post-multiplied 
matrix in such a way that the {N{m — 1) + n)th row of the post-multiplied matrix becomes 
the (2M{n — 1) + m)th row for m = 1, 2, ■ ■ ■ , 2M and n = 1, 2, ■ ■ ■ , A^, i.e., the rows having 
indexes {N{m — 1) + n), for m = 1, 2, ■ ■ ■ , 2M, are grouped for each n. Then, ^2A/,Ar^-E'2Af at 
is a (2MiV)-by-(2MiV) block diagonal matrix with block size (2M)-by-(2M), because Q, is 
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a {2MN)-hy-{2MN) block matrix with diagonal blocks of block size A^-by-iV. Similarly, the 
{MN)-hy-{MN) matrix E2m,n/2 is defined. Let Ei and E2 be the {MN)-hy-{MN) matrix that 
are defined as Ei = Pmn{In/2 ® E2m)E2m,n/2 and E2 = {In/2 ® E2m)E2m,n/2, respectively, 
where the (2M)-by-(2M) matrix E2M permutes the rows of the post-multiplied matrix in such 
a way that the (2m)th row becomes the mth row and that the (2m — l)th row becomes the 
(M + m)th row for m = 1, 2, ■ ■ ■ , M. Then, the definition S = Gm,n^G]^in in (l30al) leads to 



E2M,N^E2M^N 



EiITE^ Omn 



Omn E2TiE2 



(45) 



Since E2m,n, Ei, and E2 are all permutation matrices, e^p^ in (|44a| ) can be rewritten as il^^^ = 
tr{S*-S*(a2/2j,^Ar + S*)-iS* + S-S(fr2/2j^^^ + S)-iS}/(2MiV), which leads to ii^^ = ef^ 
because S is Hermitian symmetric. Therefore, the conclusion follows. □ 

The proposed estimator requires the pre-processing of r to obtain s, the computation and the 
multiplication of Fp, and the inverse operation on yp to obtain Xp. As it can be seen from (fT9l) and 
(|20l) . the major burden in the asymptotic FRESH properization and its inverse operation comes 
from the multiplications of the (M N)-hy-(M N) centered DFT matrix. This can be efficiently 
implemented with only C(MA^log(MA^)) complex-valued scalar multiplications by using, e.g., 
the fast algorithm in [271. As it can be seen from d39] ). the major burden in computing Fp comes 
from the inversion of {(j'^I2mn + S). The inversion of an {M N)-hy-{M N) block matrix with 
diagonal blocks of size A^-by-A^ requires 0{M^N) complex- valued scalar multiplications [j24|. 
Thus, this inversion of (cr^/2A/Af + 5]) has the same order 0{M^N) of computational complexity 
because S is an {MN)-by-{MN) block matrix with diagonal blocks of size (iV/2)-by-(iV/2). 

Since we consider the block processing where M is a fixed small number and A^ is much larger 
than M, the overall complexities of the WLMMSE and the proposed estimators as functions only 
of A^ can be rewritten now as 0{N^) and C(A^logA^), respectively. Thus, the computational 
complexity of the proposed estimator is much lower than that of the WLMMSE estimator. 
Though not major, an additional complexity reduction comes from the fact that the proposed 
estimator is linear that requires the computation and multiplication of one matrix Fp instead of 
two matrices -Fi,opt and -F2,opt, all with MN(^ 1) columns. 
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B. Numerical Results 

In this subsection, numerical results are provided that show the computational efficiency and 
the asymptotic optimality of the proposed estimator. 

The first result is to compare the complexity of the proposed estimator with that of the 
WLMMSE estimator. Fig. [9] shows that the number of complex-valued multiplications needed 
in computing the WLMMSE and the proposed estimators for cycle periods M = 1, 2, and 4. 
Recall that the computational complexities of the WLMMSE and the proposed estimators are 
0{N^) and 0{N log N) for a fixed integer M, respectively. It can be seen from Fig. |9]that the 
computational complexity of the proposed estimator is much lower than that of the WLMMSE 
estimator As predicted, the approximation motivated by Theorem |2] leads to this significant 
complexity reduction. 

The second result is to show the asymptotic optimality of the proposed estimator We consider 
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the case where an improper-complex SOCS random process is obtained by uniformly sampling 
a CT zero-mean improper-complex SOCS random process. The CT random process is, e.g., 
a Gaussian jamming signal, generated by OQPSK modulating two independent real-valued 
independent and identically distributed zero-mean symbol sequences with the SRRC pulse having 
roll-off factor 0.22. This signal is sampled at 2-times the symbol rate of the OQPSK symbols, 
which results in the DT zero-mean improper-complex SOCS random process with cycle period 
M = 2. Thus, the entries of x in (|3TI) are the MN consecutive samples of the random process. 
Fig. [To] shows that the average MSEs of the WLMMSE and the proposed estimators versus N 
for symbol energy per noise density Es/Nq = 0, 5, and 10 [dB], where the theoretical results 
are evaluated by (1361 ) and (|40|) while the simulated results are obtained from 10^ Monte-Carlo 
runs. It can be seen that, as shown in Theorem [3] the average MSE of the proposed estimator 
approaches that of the WLMMSE estimator as N increases. 
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V. Application of Asymptotic FRESH Properizer to Signal Presence Detection 

Problem 

In this section, again given a finite number of consecutive samples, now the signal presence 
presence detection of a zero-mean improper-complex SOCS Gaussian random process is consid- 
ered in additive proper-complex white Gaussian noise. It is well known ll28ll that the likelihood 
ratio is a sufficient statistic for all binary hypothesis tests under any optimality criterion. In [|5l, 
it is shown that the exact LRT statistic can be written as a quadratic function of the augmented 
observation vector when the improper-complex random process is Gaussian. Similar to the 
estimation problem in the previous section, this optimality of the use of the augmented vector 
still holds even if the output of the asymptotic FRESH properizer as an equivalent observation 
vector is processed. Again motivated by Theorem |2] we propose to use this equivalent observation 
vector of half the length of the augmented vector as if the signal component is a proper-complex 
random vector with its frequency-domain covariance matrix being equal to the masked version 
of the exact frequency-domain covariance matrix. It turns out that this suboptimal test statistic is 
asymptotically optimal in the sense that its difference from the exact LRT statistic having higher 
computational complexity converges to zero w.p. 1 as the number of samples tends to infinity. 

A. Asymptotically Optimal Low-Complexity Detector 
Let r be the length- MA^ observation vector modeled by 

Ho : r = V 

versus (46) 

"Hi : r = x + V, 

where T-Lq and Tii are the null and the alternative hypotheses, respectively. The desired signal 
X to be detected consists of the MN consecutive samples of a zero-mean improper-complex 
SOCS Gaussian random process with cycle period M eN, and v is the additive proper-complex 
white Gaussian noise vector. 

Throughout this section, the definitions and the notations of the second-order statistics of x, 
V, and r follow those in the previous section except that under the null hypothesis the second- 
order statistic of r does not contain the desired signal component. Similar to (|2T)) . the augmented 
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observation vector is denoted by f. Then, the LRT statistic of the observation vector r is provided 
as follows. 

Lemma 10: Given the observation model (|46] ). the LRT statistic of r is given by 

t/rJ = -l-f^i?sf, (47) 

LRT 2MN 

where R^ = (x~^I2mn - {(^^hhiN + ^a-)"^ 

Proof: It is straightforward by using the probability density function (PDF) of the augmented 

improper-complex Gaussian random vector [591. For details, see (51 Section 7.4]. □ 

Note that, since the length-MA^ observation vector r is improper and Gaussian, the LRT 
statistic Tj^YiT ^^ (l47l) is a quadratic form of the length-2MA^ augmented observation vector r. 
Under any optimality criterion such as the Neyman-Pearson, the Bayesian, and the minimax 
criteria, the optimal detector computes T^^^^, compares it with an optimal threshold 7], and then 
declares Hi if T^^j. > rj and Hq otherwise. 

The computation of T^^^ requires the matrix inversion of {<j'^I2mn + Rx)^ which is the 
(2MA^)-by-(2MA^) covariance matrix of the augmented observation vector r under "Hi, and 
the matrix-vector multiplication. The major burden in computing T-l^\. comes from the matrix 
inversion and it requires 0{M^N'^) complex- valued scalar multiplications. The matrix-vector 
multiplication requires only 0{M'^N'^) complex- valued scalar multiplications. Thus, the overall 
computational complexity of the LRT statistic is 0{M'^N'^). 

An altemative observation model that is equivalent to the original one in (1461 ) can be obtained 
by applying the asymptotic FRESH properization to r as 

"Ho : s = w 

versus (48) 

"Hi : s = y + w, 

where s = f{r), y = f{x), and w = f{v). Since the equivalent observation vector s is 
not proper in general, the augmentation of s and s* is still needed to compute the exact LRT 
statistic of s. Instead, we propose to use a suboptimal LRT statistic as follows, which is a 
quadratic function only of s. 

Definition 15: Given the observation model (l48T l. the proposed test statistic of s is defined by 

T^ = -^s^Wf^^llWMNS, (49) 
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where S = g^'^Imn - {ct'^Imn + S)^\ 

It can be immediately seen that Tp ' can be viewed as the exact LRT statistic of s when 



the covariance and the complementary co variance matrices of y are set equal to W/}^SWj 



MN 

and Omn appearing in (|29al ) and (|29b| ). respectively. Similar to the estimation problem in the 
previous section, this idea of using the asymptotically equivalent matrices is again motivated by 
Theorem |2l which naturally leads to the asymptotic equivalence of the LRT statistic Tl^^I in 
(1471) and the proposed test statistic T^^^ in (l49l) . 

Now, the asymptotic optimality of the proposed test statistic is provided. 

Theorem 4: The proposed test statistic approaches the exact LRT statistic as the number of 
samples tends to infinity in the sense that 

lim(Tf)-TS) = 0,w.p. 1, (50) 

under "Hq and Hi. 

Proof: We define Ct = a~'^l2MN — (c"^-f2A/Jv + ^)^ and introduce a suboptimal test statistic 

Tili ^ ^'-'"'^MN^^MNr, (51) 

which is obtained by replacing Rx in (l47l) with WIJj^CIWmn and by using the fact that Wajn is 
unitary. This introduction of T^^j, is motivated by W^/jyi?^ W^/jv ~ ^ in (|26al) . Similar to the 
proof of Theorem |3l it can be shown that Rg. ~ Wjjjy^ClWMN- In |[30| . the LRT problems 
are considered where the covariance matrix of a proper-complex Gaussian signal vector is 
either Toeplitz or block Toeplitz. It is shown under both T-Lq and Tii that a suboptimal test 
statistic, which is obtained by replacing the covariance matrix in the quadratic form of the LRT 
statistic with its asymptotically equivalent one, converges to the exact LRT statistic w.p. 1 as 
the number of samples tends to infinity [30, Propositions 1 and 3]. By applying this result, we 
have limAr^oo(^LRT ~ ^lrt) = ^-P- ^ because R^ ~ WIjj^CIWmn- Therefore, in order to 
show (|50l) . it now suffices to show T^yit = Tp . 

Let s be the frequency-domain equivalent observation vector that is defined as s = Wmns. 
Then, by using the permutation matrices E2m,n, Ei, and Ei that are defined in the proof of 
Theorem |3] in Section IV, we have 



E2m,nWmnT 



Eis* 

E2S 



(52) 
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By using ^ and ^, we can rewrite fj^^j, defined in dSB as f^^l = (s'^E*s*+s^Ss)/(2MiV) 
Tp ' because s^Ss is real-valued and (|45T ) leads 



-E'2M,7V^-E'2A'/,Ar 



EiYl*E^ Omn 

Omn E2^E2 



(53) 



Therefore, the conclusion follows. □ 

The computation of the proposed test statistic Tp requires the pre-processing of r to obtain 
s, the matrix inversion of (a'^lMN + S), and the matrix-vector multiplication to compute the 
quadratic form. As shown in Section IV, the asymptotic FRESH properization of the length- MA^ 
vector r and the matrix inversion of {(t^Imn + S) require 0{MN\og{MN)) and 0{M^N) 
complex-valued scalar multiplications, respectively. In addition, the matrix- vector multiplication 
requires 0{M'^N) complex- valued scalar multiplications, because (a'^I^N + S) is an (MN)- 
by-(MA^) block matrix with diagonal blocks of size (A^/2)-by-(A^/2). Since we consider the 
block processing where M is a fixed small number and N is much larger than M, the overall 
complexity in computing the exact LRT statistic and the proposed test statistic as functions only 
of A^ can be rewritten now as 0{N^) and C(A^logA^), respectively. Thus, the computational 
complexity of the proposed test statistic is much lower than that of the exact LRT statistic. 

B. Numerical Results 

In this subsection, numerical results are provided that show only the asymptotic optimality of 
the proposed test statistic because its computational efficiency can be similarly shown as Fig. |9l 
Throughout this subsection, the zero-mean improper-complex SOCS Gaussian random process 
to be detected is obtained by uniformly sampling the OQPSK-like signal generated in a similar 
way to that described in Section IIV-BI 

The first result is to show the convergence of the proposed test statistic to the exact LRT 
statistic. Fig. \TT\ shows the statistical average of the exact LRT and that of the proposed test 
statistics for symbol energy per noise density Es/Nq = 0, —5, and —10 [dB]. It can be seen that 
the statistical averages of the exact LRT and the proposed test statistics approach the same non- 
zero value as the number of samples tends to infinity under both hypotheses. In particular, they 
coincide under Hi because the trace of S are the same as that of R^ . Since the convergence 
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Fig. 11. Statistical averages of the exact LRT and the proposed test statistics versus A'^ for Es/Nq = 0, —5, and —10 [dB], 
under (a) 'Ho and (b) T-Li. 



w.p. 1 is hard to show by using Monte-Carlo simulations, we instead show the convergence in 
probability that is implied by the convergence w.p. 1 BUI . Fig. [T2l shows the empirical cumulative 
distribution function (CDF) of the difference of the proposed test statistic from the exact LRT 
statistic for symbol energy per noise density Es/Nq = —5 [dB] and N = 100, 200, and 400, 
where the empirical CDF is obtained by 10^ Monte-Carlo runs. It can be seen that, as implied 
by Theorem |4l the empirical CDF of the difference quickly converges to the unit step function 
as the number of samples increases under both hypotheses. 

The second result is to compare the receiver operating characteristic (ROC) curves of the 
optimal detector that uses the exact LRT statistic (l47l) and the proposed detector that uses the 
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proposed test statistic (I491 ). As a common practice in computing the PDF of a quadratic function 
of Gaussian random vectors [|32l . we approximate the statistics by gamma random variables 
to obtain the probability of miss Pm — Pr(declare "Hol^i) and the probability of false alarm 
Pfa — Pr (declare l-Li\l-Lo)- The CDF Fx{x;a,b) of the gamma random variable X with two 
parameters a and b is given by 



Fx{x;a,b) = I i—,b-l\ , Vx > 0, 
where the Pearson's form of incomplete gamma function I{u,p) is defined by 

1 fu^/p+T 



I{u,p) 



rip+l)Jo 



fPe-'dt 



(54) 



(55) 
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and where the gamma function r(x) is defined by r(x) = J^t^~^e~^dt OTIl . Note that the 
mean and the variance of the random variable whose CDF is Fx{x; a, h) are given by ah and 
a6^, respectively. Thus, Pm and Pfa of the optimal and the proposed detectors can be computed 
by using the conditional means and variances of the statistics under 1-Lq and 1-Li. Fig. [13] shows 
the ROC curves of the optimal and the proposed detectors for symbol energy per noise density 
Es/Nq = —5 [dB] and A^ = 10, 100, and 250, where the theoretical results are evaluated by using 
(|54|) while the simulated results are obtained from 10^ Monte-Carlo runs. It can be seen that 
the performance of the detectors is accurately approximated by using the gamma distributions. 
It can be also seen that the proposed detector performs almost the same as the optimal detector 
does even when A^ = 10. 



40 SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY (VER 1.0) 

VI. Conclusions 

In this paper, the asymptotic FRESH properizer is proposed as a pre-processor for the block 
processing of a finite number of consecutive samples of a DT improper-complex SOCS random 
process. It turns out that the output of this pre-processor can be well approximated by a proper- 
complex random vector that has a highly structured frequency-domain covariance matrix for 
sufficiently large block size. The asymptotic propriety of the output allows the direct application 
with negligible performance degradation of the conventional signal processing techniques and 
algorithms dedicated to the block processing of proper-complex random vectors. Moreover, the 
highly- structured frequency-domain covariance matrix of the output facilitates the development 
of low-complexity post-processors. By solving the signal estimation and signal presence detection 
problems, it is demonstrated that the asymptotic FRESH properizer leads to the simultaneous 
achievement of computational efficiency and asymptotic optimality. Further research is warranted 
to apply this pre-processor to various communications and signal processing problems involving 
the block processing of a DT improper-complex SOCS random process and to derive such almost 
optimal low-complexity post-processors. 

References 

[1] F. D. Neeser and J. L. Massey, "Proper complex random processes with applications to information theory," IEEE Trans. 

Inf. Theory, vol. 39, no. 4, pp. 1293-1302, July 1993. 
[2] P. J. Schreier and L. L. Scharf, "Second-order analysis of improper complex random vectors and processes," IEEE Trans. 

Signal Process., vol. 51, no. 3, pp. 714-725, Mar. 2003. 
[3] P. Wahlberg and P. J. Schreier, "Spectral relations for multidimentional complex improper stationary and (almost) 

cyclostationary processes," IEEE Trans. Inf. Theory, vol. 54, no. 4, pp. 1670-1682, Apr. 2008. 
[4] G. Taubock, "Complex-valued random vectors and channels: Entropy, divergence, and capacity," IEEE Trans. Inf. Theory, 

vol. 58, no. 5, pp. 2729-2744, May 2012. 
[5] P. J. Schreier and L. L. Scharf, Statistical Signal Processing of Complex-Valued Data: The Theory of Improper and 

Noncircular Signals. NY: Cambridge Univ. Press, 2010. 
[6] B. Picinbono and P. Chevalier, "Widely linear estimation with complex data," IEEE Trans. Signal Process., vol. 43, no. 8, 

pp. 2030-2033, Aug. 1995. 
[7] W. M. Brown and R. B. Crane, "Conjugate linear filtering," IEEE Trans. Inf. Theory, vol. 15, no. 4, pp. 462-465, July 

1969. 
[8] E. Serpedin, F. Panduru, I. Sari, and G. B. Giannakis, "Bibliography on cyclostationarity," Signal Process., vol. 85, no. 

12, pp. 2233-2303, Dec. 2005. 
[9] W. A. Gardner, A. Napolitano, and L. Paura, "Cyclostationarity: Half a century of research," Signal Process., vol. 86, no. 

4, pp. 639-697, Apr. 2006. 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY (VER 1.0) 41 

[10] W. A. Gardner and L. E. Franks, "Characterization of cyclostationary random signal processes," IEEE Trans. Inf. Theory, 

vol. 21, no. 1, pp. 4-14, Jan. 1975. 
[11] W. A. Gardner, "Signal interception: A unifying theoretical framework for feature detection," IEEE Trans. Commun., vol. 

36, no. 8, pp. 897-906, Aug. 1988. 
[12] W. A. Gardner and C. M. Spooner, "Signal interception: Performance advantages of cyclic-feature detectors," IEEE Trans. 

Commun., vol. 40, no. 1, pp. 149-159, Jan. 1992. 
[13] W. A. Gardner, W. A. Brown, and C.-K. Chen, "Spectral correlation of modulated signals: Part II-Digital modulation," 

IEEE Trans. Commun., vol. 35, no. 6, pp. 595-601, June 1987. 
[14] W. A. Gardner and W. A. Brown, "Frequency- shift filtering theory for adaptive co-channel interference removal," in Proc. 

23rd Asilomar Conf. Signals, Systems, Computers, vol. 2, Pacific Grove, CA, 30 Oct.- 1 Nov. 1989, pp. 562-567. 
[15] G. B. Giannakis, "Filterbanks for blind channel identification and equalization," IEEE Signal Process. Lett., vol. 4, no. 6, 

pp. 184-187, June 1997. 
[16] A. Sabharwal, U. Mitra, and R. Moses, "MMSE receivers for multirate DS-CDMA systems," IEEE Trans. Commun., vol. 

49, no. 12, pp. 2184-2197, Dec. 2001. 
[17] J. H. Cho, "Joint transmitter and receiver optimization in additive cyclostationary noise," IEEE Trans. Inf. Theory, vol. 50, 

no. 12, pp. 3396-3405, Dec. 2004. 
[18] Y. H. Yun and J. H. Cho, "An optimal orthogonal overlay for a cyclostationary legacy signal," IEEE Trans. Commun., vol. 

58, no. 5, pp. 1557-1567, May 2010. 
[19] B. W. Han and J. H. Cho, "Capacity of second-order cyclostationary complex Gaussian noise channels," IEEE Trans. 

Commun., vol. 60, no. 1, pp. 89-100, Jan. 2012. 
[20] W. A. Gardner, "Cyclic Wiener filtering: Theory and method," IEEE Trans. Commun., vol. 41, no. 1, pp. 151-163, Jan. 

1993. 
[21] B. W. Han and J. H. Cho, "Optimal signaling in second-order cyclostationary Gaussian jamming environment," in Proc. 

of IEEE Military Communications Conf (MILCOM), Baltimore, MD, 7-10 Nov. 2011, pp. 102-107. 
[22] J. Eriksson and V. Koivunen, "Complex random vectors and ICA models: Identifiability, uniqueness, and separability," 

IEEE Trans. Inf Theory, vol. 52, no. 3, pp. 1017-1029, Mar. 2006. 
[23] R. M. Gray, "Toeplitz and circulant matrices: A review," Found. Trends Commun. Inf. Theory, vol. 2, no. 3, pp. 155-239, 

2006. 
[24] Y G. Yoo and J. H. Cho, "Asymptotically optimal low-complexity SC-FDE in data-like co-channel interference," IEEE 

Trans. Commun., vol. 58, no. 6, pp. 1718-1728, June 2010. 
[25] J. H. Yeo and J. H. Cho, "Properization of second-Order cyclostationary random processes and its application to signal 

presence detection," in Proc. IEEE Military Communications Conf. (MILCOM), Orlando, FL, 29 Oct.-l Nov. 2012, pp. 

1-6. 
[26] B. Picinbono and P. Bondon, "Second-order statistics of complex signals," IEEE Trans. Signal Process., vol. 45, no. 2, 

pp. 411-420, Feb. 1997. 
[27] G. Bi and Y Chen, "Fast generalized DFT and DHT algorithms," Signal Process., vol. 65, no. 3, pp. 383-390, Mar. 1998. 
[28] H. V. Poor, An Introduction to Signal Detection and Estimation, 2nd ed. NY: Springer- Verlag, 1994. 
[29] A. van den Bos, "The multivariate complex normal distribution-A generalization," IEEE Trans. Inf. Theory, vol. 41, no. 

2, pp. 537-539, Mar. 1995. 



42 SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY (VER 1.0) 

[30] W. Zhang, H. V. Poor, and Z. Quan, "Frequency-domain correlation: An asymptotically optimum approximation of quadratic 
likelihood ratio detectors," IEEE Trans. Signal Process., vol. 58, no. 3, pp. 969-979, Mar. 2010. 

[31] A. Papoulis, Probability, Random Variables, and Stochastic Processes, 4th ed. New York: McGraw-Hill, 2002. 

[32] U. Grenander, H. O. Pollak, and D. Slepian, "The distribution of quadratic forms in normal variates: A small sample 
theory with applications to spectral analysis," J. Soc. hid. Appl. Math., vol. 7, no. 4, pp. 374-401, Dec. 1959. 



